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Colin Benjamin and Jiannis K. Pachos 
Quantum Information Group, School of Physics and Astronomy, 
University of Leeds, Woodhouse Lane, Leeds LS2 9JT, UK. 

We propose to combine the advantages of graphene, such as easy tunability and long coherence 
times, with Josephson physics to manufacture qubits. If these qubits are built around a and n 
junction they can be controlled by an external flux. Alternatively, a d-wave Josephson junction can 
itself be tuned via a gate voltage to create superpositions between macroscopically degenerate states. 
We show that ferromagnets are not required for realizing n junction in graphene, thus considerably 
simplifying its physical implementation. We demonstrate how one qubit gates, such as arbitrary 
phase rotations and the exchange gate, can be implemented. 



I. INTRODUCTION 

Graphene, a monatomic layer of graphite exhibits 
promising electronic properties that can be employed for 
quantum technologies 1 . Characteristically, its low energy 
excitations are described by the Dirac equation, it has a 
zero band gap, electronic speeds can reach a hundredth 
of the speed of light and it supports long range phase 
coherence. However it has not yet been utilised to create 
qubits suitable for quantum computation, apart from a 
proposal which meshes it with bilayer structures^. Here 
we show that a key ingredient of Josephson qubits, a 
7T— junction^ can be easily generated in graphene by ap- 
plication of a gate voltage alone. We establish a para- 
metric regime for observing this effect and show how to 
manufacture qubits. These Josephson qubits can be used 
to perform single quantum gates, such as the phase and 
exchange gates. This opens up the possibility of employ- 
ing graphene and utilizing its advantages for quantum 
information processing^. 

The physical system we employ consists of a graphene 
substrate with superconducting correlations induced in 
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sections via the proximity effect^ or by turning graphene 
superconducting^ via doping. It comprises of two d- 
wave Josephson junctions (distinguished by their ground 
states, one at a phase difference </> = and the other at 
cf) = 7r), arranged as in Fig. 1(a). The total energy of the 
system is controlled by the flux, <!>, that passes through 
the ring. The reversal of super-current in a Josephson 
device, where the free energy has global minima at phase 
difference <f> = n, is referred to as ir shift. The corre- 
sponding Josephson junction is termed a tt junction. This 
is in contrast to a junction wherein the free energy has 
a global minimum at phase difference <ft = 0£. To be able 
to encode a qubit we have to construct a ir junction and 
integrate it with the rest of our device (the junction). 
A 7r junction is needed to create a doubly degenerate 
ground state, where a qubit is encoded. Here we demon- 
strate that a 7T junction can be identified in our system 
without the need of any ferromagnetic elements^, thus 
greatly simplifying its experimental implementation. In 
Fig. 1(b) we depict a simple d-wave graphene Joseph- 
son junction, which has two degenerate ground states, 
that can encode a qubit. In particular, we prove that a 
complete set of single qubit gates can be efficiently im- 
plemented demonstrating that our proposal is promising 
for quantum computation. 

In Fig. [2j we show our graphene ir junction set-up. It 
is known that with s-wave superconductors a it junction 
is not possible^. However, a Josephson junction with 
d-wave superconductors can exhibit a ir shifts. Thus, 
we consider d-wave correlations in the superconducting 
segments (see Fig. I). 



II. THEORY 



FIG. 1: (Color online) An overview of the set-up. (a) Two 
semicircular d-wave superconducting graphene strips (Gs) 
with normal graphene layers on top and bottom enclosing 
a magnetic flux, $. By the application of suitable gate volt- 
ages to the normal graphene strip the junctions are tuned 
to either <f> = or 7r phase shift, (b) A graphene d-wave 
Josephson junction. For relatively small intervening length 
between the superconducting graphene one can have situa- 
tions wherein degenerate ground states are formed and are 
pliable to external control via a gate voltage. 



The kinematics of quasi-particles in graphene is de- 
scribed by the Dirac-Bogoliubov-de Gennes equation^, 
which assumes the form 

(H-Eri AI )v = e* (i) 

where E is the excitation energy, A is the superconduct- 
ing gap of a d-wave superconductor, 'J is the wavefunc- 
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FIG. 2: (Color online) The Furusaki-Tsukada approach and 
the processes involved. Top: 8g is the angle of incidence of 
electron- like quasiparticle, while —0% is the angle of its reflec- 
tion. Hole-like quasiparticle are Andreev reflected at angle 
6$. In the normal region electron and holes are transmitted 
and incident with angles 6 and 8 A . Bottom: In type 1 process 
an electron-like quasiparticle is incident from the left, while 
in type 2 process a hole-like quasiparticle in incident from 
the left, oi, 62, d\ and d,2 are amplitudes of hole-like quasi- 
particle, while d2, 61, Ci and C2 are scattering amplitudes for 
electron- like quasi-particles. 



tion and • represents 4x4 matrices. In the above equation 

H = ( H + H± = -inv F {<T x d x ±a y d y ) + U. (2) 

Here h,v F (set equal to unity hence forth) are the 
Planck's constant and the energy independent Fermi ve- 
locity for graphene, while the cr's denote Pauli matrices 
that operate on the sub-lattices A and B. The electro- 
static potential U can be adjusted independently via a 
gate voltage or doping. We assume U = 0, in the normal 
region, while U = — Uq in the superconducting graphene. 
In our work we consider U$ = 100A. Further we choose 
d-wave superconducting correlations which imply a type 
II (or high T c ) superconductor. This is most likely to be 
observed in grapheneii. The subscripts of Hamiltonian 
± refer to the Fermi points K + and K_ in the Brillouin 
zone. T = —r y ® a y C, (C being complex conjugation) is 
the time reversal operator, with r being Pauli matrices 
that operate on the ± space and I is the identity matrix. 

To calculate the Josephson supercurrent, Free energy 
and show the formation of a 1 junction we proceed 
by first calculating the scattering wave functions of 
our system. Let us consider (Type 1 scenario in 
Fig. 2)an incident electron-like quasiparticloiS. from 
the left superconductor with pairing gap A(8 + )e 1 ^ 
(x < 0) and energy E. For a right moving electron-like 
quasiparticle with an incident angle 9 the eigen- 
vector and corresponding momentum read = 

[u e , u e e ie+ , Vee-tt , v e e i ( e+ -'t l t)] T e i i e cos e+x , q e 

(E F + U + \/E 2 - |A(6>+)| 2 ). A left moving 



electron-like quasiparticle is described by the sub- 
stitution 9 — > 7T — 9. If Andreev-reflection takes 
place, a left moving hole-like quasiparticle is gen- 
erated with energy E, angle of reflection 6~ and 
its corresponding wavefunction is given by ^ r g 1 _ = 

[v h , -v h e~ i6 ~ , u h e-^ , -utfi-< e ~ + ^\ T e- iqh 
(E F +U - yjE 2 - |A(6»-)| 2 ). The quasi-particle wave- 
vectors can also be expressed as q e l h = Ep + Uq ± l/£, 
where £ is the coherence length. For the Dirac- 
Bogoliubov dc Gennes equations to hold the Fermi 
wavelength in the superconductor 1/(E F + Uo) should 
be much smaller than the coherence length. The 
superscript e (h) denotes an electron-like (hole-like) 
excitation. Since translational invariance in the 
y-dircction holds the corresponding component of 
momentum is conserved. This condition allows for 
the determination of the Andreev reflection angle 9~ 
through q h sm(6~) = q e sm(9 + ). The coherence factors 

are given by u e/h 

defined 9 + = 9 e Sl 9~ = n — 6$, where the angles 



(l + Vl-|A(0±)|7i? 2 )/2, 
We have also 



arc defined in Fig. [5] In our study we have d-wave 
superconductors, thus A(0 ± ) = Acos(20 ± — 27) and the 

macroscopic phase is e 1 ^ 1 / 2 = e 1 ^ 1 ^ 2 



|A(e±)| 



We choose 



the superconductor oriented along the 110 direction, 
implying 7 = 7r/4. 

In the normal region the eigenvector and correspond- 
ing momentum of a right moving electron with an in- 
cident angle 9 read: ^+ = [1, e l6 , 0, Q] T e ip ° cob6x , p e = 
(E + E F ). A left moving electron is described by the sub- 
stitution 9 — > 7r — 9. If Andreev-reflection takes place, a 
left moving hole is generated with energy E, angle of re- 
flection 9a and its corresponding wave function is given 
by- 1/^ = [0,0,l,e- ie ' 4 ] T e- i P" cose - 4;E , p h = (E - E F ). 
The transmission angles 9 and 9a for the electron-like and 
hole-like quasi-particles are given by q e sin 9 e s = p e sin 9 
and q e sin 9 e s = p h sin 9a- 

The full wave function in the type 1 scenario can be 
written as below for the various regions 
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(3) 



Matching the wave functions at the interfaces one can 
solve for the amplitudes of reflection 01, i>i, c\ and d\. 
Similarly, one can write the wave functions in case of type 
2 scenario (hole incident from the right) and calculate the 
amplitudes a<i, 62, C2, and c?2- The detailed balance for 
the amplitudes are verified as follows 

Ca^E) = C'a 2 (-4>,E), 

bi(4>,E) = bi(-(t>,E)(i = 1,2), (4) 

with C = |A(e'-)| cos^si an( i C = A("e'+)| cos $s- Follow- 
ing the procedure established in RefJ^ and employing 
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analytic continuation E — > iw n the dc Joscphson current 
is calculated as 

n ra+ 



m - ZmL 1 



ax(9 + ,<t>, iw n ) 
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-]cos((9f)d(9|, 
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(5) 



where /3 = l/fc B T, fi n ,± = y/wl + |A((9±)| 2 and w n = 
Trk B T{2n + 1), n = 0, ±1, ±2, .... 

The above equation has a simple physical interpreta- 
tion^. Andreev reflection is equivalent to the breaking 
up or creation of a Cooper pair. The scattering amplitude 
a\ describes the process in which an electron-like quasi- 
particle coming from the left superconducting graphene 
strip (x < 0) is reflected as a hole-like quasiparticle. The 
amplitude a 2 corresponds to the reverse process in which 
a hole-like quasiparticle is reflected as an electron-like 
quasiparticle. This implies that a\ and a 2 correspond to 
the passage of a Cooper pair to the left and right respec- 
tively hence, the dc Josephson current is proportional 
to a\ — 02- Further, the dc Josephson current is an odd 
function of the phase difference, (j>, as seen by the detailed 
balance condition, a 2 {4>, iw n )/C = Oi(— 4>, iw n )/C. To 
calculate the Josephson current one thus takes the differ- 
ence between the amplitudes a\ and a 2 and then sums 
over the energies. In this approach we account for all the 
energies both bound states and the continuum. Eq[5]can 
be simplified as- 
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(6) 



In EqG3 A,B,A',B', and C are functions of 
9 + ,iw n , Ef&nd d. The Free energy of the Josephson 
junction can then be calculated as 

m = ^- 1* ww. (7) 

27T Jo 



III. 7r- JUNCTION 

Now we illustrate the results for the Josephson cur- 
rent as function of the length of the normal graphene 
interlude as well as the phase difference across the two 
superconducting graphene strips. The calculations are 
performed by treating Eqs. ([5]) and ([7]) numerically and 
the derived results hold for the T — ► temperature limit. 
Fig. EJa) shows the Josephson current as function of the 
Fermi energy, in the normal graphene strip, for different 



lengths of the normal graphene layer. Note that Fermi 
energy is easily controllable in graphene. The plot shows 
that for extremely small length of normal graphene layer 
the Josephson current is negative for a wide range of 
Fermi energy, implying a tt shift, while for larger in- 
tervening normal layers the Josephson current changes 
sign at larger values of the Fermi energy One important 
fact to note is that for increased d the current decreases, 
which is in agreement with past Joscphson works. An- 
other observation from Fig. G2a) is that at large Fermi 
energy the Josephson supercurrent becomes independent 
of Ef. The explanation for this is- when Ep >> E,A, 
the angles for electron and hole-like quasi-particles are 
9g = 9g = 9 = —9 a- With this condition, the factor J 
from Eq|6j the Josephson supercurrent shorn of all pref- 
actors, reduces to- 
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In the above equation, 7 = (p e + ph)dcos(9) = 
Edcos(9),h = (E - x)/2E,g = (E + x)/2E,x = 

E 2 — sin(6>) 2 . Thus in this limit the Josephson su- 
percurrent becomes completely independent of Ep. Fur- 
ther, for d — > one can clearly see from Fig. [3](a) that 
the Josephson supercurrent becomes completely nega- 
tive, this is also evident from Eqj8l wherein J reduces 
to -2w„sin(26»)/(2u; I 2 l + sin(29) 2 ),E = iw n . Fig. G^b) 
shows the current-phase relation for two different values 
of the Fermi energy. It again confirms the earlier indica- 
tion of 7r shift. Finally, to establish beyond doubt that 
as function of Fermi energy one generates a it junction 
we plot the free energy in Fig. [3fc). The plot shows that 
as one changes the Fermi energy via a gate voltage one 
changes the ground state of the junction from to it. 

As shown in Fig.[3^c-d), the Free energy, F, has a min- 
imum at (/> = 7r (for the 7r junction case) and the variation 
of F with (j> is strongly dependent on the length d and the 
Fermi energy. In this parameter regime the Free energy 
can be approximated as F ~ — £ l 7r [ cos (07r + t) + 1], with 
E^ being the Josephson coupling constant. The and 
7r junctions, depicted in Fig. 1, have Josephson energies 
U = E \ sm(</>o/2)| and U v = —E^[co8{(j)^ +7r) + l] plot- 
ted in Fig. [3fd) . The superconducting phase difference 
is <fio for the junction and 0^ for the ir junction. The 
total flux in the ring <i> satisfies (f) n — <f)Q = 27r$/$ ~ 27rZ, 
where $0 is the flux quantum and I is an integer. 



IV. QUBITS AND GATES 

In Ref. the authors demonstrate a qubit with a tt 
(SFS) junction^ and a (SNS) junction coupled into a 
ring. In our work we predict that our graphene based 
system, which does not need any ferromagnetic element 
in contrast to ReffTl. could implement a qubit. Further 
we show how to implement single qubit gates using our 
set up. The full Hamiltonian of the graphene ring system 
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FIG. 3: (Color online) (a) Current (in units of eA/fi and nor- 
malized by 1//3 throughout in all succeeding figures) versus 
Fermi energy, Ef, at phase difference <f> = n/2, for differ- 
ent values of width <f(in units of Hvf/A), Uo = lOO&fcsT = 
0.0001 in this and all succeeding figures, (b) Current ver- 
sus phase, where the length of normal graphene strip is 
d — 0.1, the dashed (RED) line is multiplied by a factor of 
ten for better visibility, (c) Free energy (normalized by 1//3) 
of Gs — Gn — Gs junction versus phase difference for differ- 
ent Fermi energies with junction (Ef = 2000 red dashed 
line) and n junction (Ef = fOO black solid line) and length 
of normal graphene strip d = O.f . (d) The approximate forms 
for the and ir junction energies are in good agreement with 
the real free energies and are used in analyzing the graphene 
Josephson qubit. 



(Fig. 1) is given by H = K + U to t with U tot = U + 
U„ + Ul, where Ul = — $ ext ) 2 /2Ls is the magnetic 
energy stored in the ring and K is the flux independent 
kinetic energy. We next minimize the Hamiltonian with 
respect to flux and obtain Q((f>o) — /3$o sin(0 x ) + $> ex t, 
with (3 = 2ttEt T Ls/^q- Substituting this equation in the 
expression for Utot, we have: 



tot n . 

= a\ sin — 



7T$ 
$0 



)|] + [cos(^r) - 1] +7r/3sin 2 (^). 

(9) 

with a = Eq/Ett. For typical values mentioned in Fig.|4j 
we plot Eq. ©. We observe that the energy has dou- 
ble minima located approximately at </>„■ ~ 3n/5 (|0) 
state) and 7ir/5 (|1) state) which form the basis of the 
qubit. For single layer graphene with junction area— 
0.8 x 10~ 12 /i m 2 and depth 1 nm, the electrostatic en- 
ergy E c is 2.5 x 10~ 24 J, while Eq the junction energy for 
the zero junction is around 1000i? c . Thus for a = 3.0, we 
have AE, the energy gap, between the ground and first 
excited state AE/h = 1000GHz. The basic phase gate 



FIG. 4: (Color online) (a) Normalized energy, Utot/E w , as 
function of 4>-n for no external magnetic flux, (b) In presence 
of an external magnetic field with a — 2.5. (c) The degenerate 
ground states of a d-wave Josephson junction. 



with (j) = AEAt/h = 7T could be implemented with gate 
time At given by 1 pico-second. In Fig. 4(c), the Free 
energy of a basic d-wave graphene Josephson junction is 
plotted for different values of Fermi energy and width 
d = 0.001. One can easily see that degenerate states are 
formed at </> ~ 7r/2 and 37r/2. The coupling between these 
states can be easily varied by the gate voltage effectively 
realizing single qubit gates as aforementioned. 

We will now show how to implement an exchange gate 
a x acting on the qubit states |0) and |1) for the struc- 
ture as depicted in Fig. 1(a). This is realized by a tun- 
nelling transition between the potential minima that en- 
code these qubit states. Assuming the coupling potential 
is deep enough we approximate the qubit states by Gaus- 
sians centered at the minima of Utot- By varying a (or 
E c ) one can induce tunnelling between the two minima 
in a controlled way. The exchange coupling of our system 
is calculated as 



J 



: (07r-0|O>)(- 



-AE r 



(10) 

In Fig. [5] we plot the exchange coupling versus the nor- 
malized Josephson energy for various values of the elec- 
trostatic energy, E c in units of E n . We see that for 
large a no tunnelling occurs, while for a ~ 3.0 we ob- 
tain J ~ lO^-EV (for E c = 0.01) and, thus, the a x gate 
can be implemented in At ~ 10~ 6 seconds. 

To conclude we have shown a novel implementation of 
a Josephson qubit using graphene as a substrate. Our 
work is the first to predict a qubit using only monolayer 
graphene. It was shown that a ferromagnetic graphene 
layer is unnecessary to create a 7r-shift, a completely novel 
result. 7T junctions have special role in a host of appli- 
cations ranging from their use in superconducting digi- 
tal circuits to superconducting qubits. We have shown 
how a 7T junction is formed in graphene where it can 
be very easily tuned by the application of a gate volt- 
age alone. Secondly, we propose Josephson qubits and 
we present the phase and exchange gates for quantum 
computation purposes. Future proposals to make CNOT 
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FIG. 5: (Color online) Exchange coupling J (in units of 
as function of a for different values of E c . Here E c <C E^ . 



or other two-qubit gate designs could also be envisaged 
using the above architecture. 
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